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Abstract 

When a potential for a scalar field has two local minima there arise spherical 
shell-type solutions of the classical field equations due to gravitational attraction. 
We establish such solutions numerically in a space which is asymptotically de Sitter. 
It generically arises when the energy scale characterizing the scalar field potential is 
much less than the Planck scale. It is shown that the mirror image of the shell appears 
in the other half of the Penrose diagram. The configuration is smooth everywhere 
with no physical singularity. 


1. Introduction 

Gravitational interactions, which are inherently attractive for ordinary matter, can 
produce soliton-like objects even when such things are strictly forbidden in flat space. 
They are possible as a consequence of the balance between repulsive and attractive forces. 
One such example is a monopole or dyon solution in the pure Einstein-Yang-Mills theory in 
the asymptotically anti-de Sitter space BSi- In the pure Yang-Mills theory in flat space 
there can be no static solution at all Q but once gravitational interaction is included there 
arise particle-like solutions [Q. Whereas all solutions are unstable in the asymptotically 
flat or de Sitter space, there appear a continuum of stable monopole and dyon solutions in 
the asymptotically anti-de Sitter space. The stable solutions are cosmological in nature; 
their size is typically of order where A is the cosmological constant. 

The possibility of false vacuum black holes has also been explored. Suppose that the 
potential in a scalar field theory has two minima, one corresponding to the true vacuum 
and the other to the false vacuum. If the universe is in the false vacuum, a bubble of the 
true vacuum is created by quantum tunneling which expands with accelerated velocity. 
The conhguration is called a bounce @]. Now flip the conhguration ||^. The universe 
is in the true vacuum with potential V = 0 and the inside of a sphere is excited to the 
false vacuum with V > 0. Is such a de Sitter lump in Minkowski space possible? If the 
lump is too small it would be totally unstable. The energy localized inside the lump can 
dissipate to spatial inhnity. If the lump is big enough the Schwarzschild radius becomes 
larger than the lump radius so that the lump is inside a black hole. The energy cannot 
escape to inhnity. It looks like a soliton in Minkowski space. However, as a black hole it 
is a dynamical object. The conhguration is essentially time-dependent. This false vacuum 
black hole conhguration, however, does not solve the equations at the horizon. It has been 
recently proven that there can be no such everywhere-regular black hole solution i I- 
Rather, false vacuum lumps in hat space evolve dynamically 0, MM- 

The purpose of this paper is to report new solutions to the coupled equations of grav¬ 
ity and scalar held theory which display spherically symmetric shell structure |^. We 
demonstrate that such structure appears when the potential for a scalar held has two local 
minima and the space is asymptotically de Sitter. In the examples we present, both the 
inside and outside of the shells are de Sitter space with the same cosmological constant. 
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The structure becomes possible only when the energy scale of the scalar held potential 
becomes small compared with the Planck scale. While such shell structures might not 
be easy to create in the present universe, it is quite plausible that they could have been 
created during a phase transition early in the universe 0 . A similar conhguration has 
been investigated in ref. [|H . 


The plan of our paper is as follows. In section 2 we precisely state the problem, solve 
the held equations in those regions of space-time where they can be linearized, and sketch 
the solution in the nonlinear shell region in static coordinates. In section 3 we solve the 
nonlinear equations in the shell region and display the dependence on the parameters of 
the theory. In section 4 we extend the solution from static coordinates, which have a 
coordinate singularity, to other coordinate systems that do not, thereby displaying the 
existence and character of the solution throughout the full space-time manifold. In section 
5 we study the stability of the classical solution to quantum fluctuations. Our summary 
and conclusions are given in section 6 . 


2. Shell Solutions in Static Coordinates 

Consider a scalar held coupled to gravity with the Lagrangian 


( 2 . 1 ) 


where the potential ld[0] has two minima at /i and /2 separated by a barrier. See hg. |l|. 

We look for spherically symmetric conhgurations in which the metric of space-time is 
written as 


H 

H 


+ r^{d6^ + sirr Odp ) . ( 2 . 2 ) 


The functions 0, id, and p depend only on r and t. A tetrad basis is chosen, in a region 
id > 0 , as 

Co = dt , Cl = dr , 62 = rdO , 63 = r sin 6dp . (2.3) 

P VH 

The components of the energy-momentum tensor in the tetrad basis. Tab = ^a^Gb’^Tnu, are 


2\H 
2\H 






f2 
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Figure 1: The potential V[(j)\. 


T22 = Tsj = i - V[ 4 >] 

T„i = . (2.4) 


Here dot and prime indicate t- and r-derivatives, respectively. 
The scalar held satishes 


P-. 


(9 / p • \ p d 


dt\H J r‘^dr\ p 
We introduce the integrated mass function M(t, r) by 

2GM 


0 ' + V'[(l)] = 0 . 


H = 1 


The Einstein equations are 


M = J ATir^dr Too , 

H 

— = —87iGr^(j)' , 

H 

p ( d /pH\ d / H' 2Hp' \ "I 1 — H 

2 I ^ \ p p^ / J ^ 



(2.5) 

( 2 . 6 ) 

(2.7) 

(2.8) 

(2.9) 

( 2 . 10 ) 
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One of the equations, Eq. (| 2 . 1 (J|), is redundant as it follows from Eqs. (|2.5|), (|2.7|), 


and (|2.9|) . 

We shall seek static solutions for which the set of equations reduces to 

0"(r) + reff(r)0'(r) = ^ V'[(f)] , 


M{r) = J + 1/[0]| , 


where 


eff 


= - + AirGrcj)''^ + ^ 
r H 


( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 


Eq. ( p.ll ) can be interpreted as an equation for a particle with a coordinate 0 and time r. 
Except for a factor 1/H this particle moves in a potential U[(f)] = —!/[<;/)]. The coefficient 
reff(r) represents time (r)-dependent friction. 

The potential E[0] is supposed to have two minima, at /i and / 2 . We are looking for 
a solution which starts at 0 ~ /i, moves close to 02 , and comes back to /i at r = cx). The 
particle’s potential f/[0] has two maxima. The particle begins near the top of one hill, rolls 
down into the valley and up the other hill, turns around and rolls down and then back up 
to the top of the original hill. This is impossible in flat space, as Teff is positive-dehnite so 
that the particle’s energy dissipates and it cannot climb back to its starting point. 

In the presence of gravity the situation changes. The non-vanishing energy density can 
make H a decreasing function of r so that Teff becomes negative. The energy lost by the 
particle during the initial rolling down can be regained on the return path by negative 
friction, or thrust. Indeed, this happens. 

Let us set up the problem more precisely. We take a quartic potential with /i < 0 < 02 ; 


V'[4)\ = A0(0-/i)(0-/2) , 

VW = \('t> - h){4>^ - l(k + 4/i)/ - 02(/2 - 2fi)(4, + /2)} . (2.14) 

Here V[f 2 ] = 0, and H[0] is a local maximum for the barrier separating the two minima. 
Dehne / = (|/i| -|- / 2)/2 and A/ = 02 — \ fi\- In case A/ > 0, 0 = /i corresponds to a false 
vacuum with the energy density e = V[fi] = |A/^A/ > 0, whereas 0 = 02 corresponds to 
a true vacuum with a vanishing energy density. As we shall discuss in detail below, the 
positivity of the energy density e plays an important role for the presence of shell structure. 
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but the vanishing V[f 2 ] is not essential as we see below. In a more general potential it 


could be that V[f 2 ] > V[fi]. 

We look for solutions with 0 starting at the origin r = 0 very close to /i. There is only 
one parameter to adjust: 0o = 0(0). The behavior of a solution near the origin is given by 




p = 1+ H-, 

M = + • • •, 

H = 1 — 2Gm^r‘^ + ■ 


Pa = —47rG02 , 
“3 = Y ’ 


(2.15) 


Given 0o the equations determine the behavior of a conhguration uniquely. For most 


values of 0o the corresponding conhgurations are unacceptable. As r increases, 0(r) either 


approaches 0 (the local maximum of ld[0]) or goes to —oo. Other than the two trivial 
solutions, corresponding to the false and true vacua, we have found a new type of solution. 

There are four parameters in the model, one of which, the gravitational constant 
G = rup^, sets the scale. The other three are A, /i, and 02 or, equivalently, the three 
dimensionless quantities f/mp, A///, and A. We have explored only a limited region in 
the parameter space. The moduli space of solutions depends critically on f /mp and Af /f, 
but seems to depend little on A. Nontrivial solutions appear as f /mp becomes small. 

If 00 < /i, 0(^) monotonically decreases as the radius increases to diverge to —oo. On 
the way H{r) crosses zero. Suppose instead that /i < 0o < 0 and 0o is not too close 
to /i. In the particle analogue, the particle starts to roll down the hill under the action 
of t/[0] = —F[0]. It approaches 0 = 0, and oscillates around it. In the meantime H{r) 
crosses zero. 

Now suppose that 0o is very close to, but still greater than, fp. 0(0) = /i + ^0(0) with 
0 < 50(0)// -C 1. A schematic of the resulting solution is displayed in £g. || and £g. H. 
We divide space into three regions in the static coordinates: region I (0 < r < Ri), region 
II {Ri < r < R2), and region III {R2 < r). It turns out that 0(r) varies little from fi in 
regions I and III so that the equation of motion for 0 may be linearized in those regions. In 
region II the held deviates strongly and the full set of nonlinear equations must be solved 
numerically. This is the region in which we shall hnd shell structure. H{r) deviates from 
the de Sitter value signihcantly. In region III the spacetime is approximately de Sitter 
again. H{r) crosses zero at th- 
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Figure 2: Schematic behavior of 4>{r). 


In region I the space-time is approximately de Sitter. 

...2 I T 


Too — ^ 1 hr — 1 


a = 


SirGe 


, P=l ■ 


( 2 . 16 ) 


The equation for (j){r) can be linearized with 0(r) = fi + 5(j){r). In terms of z = 


+ ( 


3 _ 5 \ d _ 

2 ~ 2 V ^ ~ 4 


/> (50 = 0 , 


(2.17) 


where = V"[fi]. This is Gauss’ hypergeometric equation. The solution which is regular 
at r = 0 is 


d0(r) = (50(0) ■ F(f + i/t, f - 1 


(2.18) 


where 


K, = 


9 

4 ’ 


2 2 

oj a = 


9mi 


Stt/A/ 'v 2 / J 


(2.19) 


We shall soon see that a solution with shell structure appears for cua 3> 1 with a particular 
choice of (50(0). The ratio of (50'(r) to (50(r) is given by 


(50'(r) 

d0(r) 


4r 


/^2 I 9 \ h’(4+m,4 _ 2r 

r ^ 1 Fr 3 ^ 3 _ 3. ^ ^ ’ • 


m, ^;z) 


( 2 . 20 ) 
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Figure 3: Schematic behavior of H{r). 


The deviation from fi at the origin, 50(0), needs to be very small for an acceptable 
solution. The behavior of the hypergeometric function for k S> 1 and 0 < z < 1 is given 
by in 


F{a + m, a — m, c; z) 


r(c) i_ _£_^i 

2; 2 + 4 


27 ? 


(1 


2;)2 4 “ exp|2Ksin ^-y/ij . (2.21) 


The ratio 50(r)/50(O) grows exponentially as r increases like — 

exp|2fi:sin“^ 7?}- of region I, 50/|/i| needs to be very small for the 

linearization to be valid. The ratio of F\z) to F{z), J{z) in (p2.20|) , is given by 

K 


J{z) = 




forKS>l, 0<z<l. 


( 2 . 22 ) 


In region II, 0(r) varies substantially and the nonlinearity of the equations plays an 
essential role. In this region the equations must be solved numerically. With hne tuning 
of the value of 50(i?i) nontrivial shell solutions will be found. 

The algorithm is the following. First 50(i?i) is chosen and 50'(i?i) is evaluated by 
( p.20| ) and ( |2.22|) . To the order in which we work the metric is H{Ri) = 1 — (i?i/a)^ and 
p{Ri) = 1. With these boundary conditions Eqs. (|2.11|) and (|2.12|) are numerically solved. 

The behavior of solutions in region II is displayed in Eg. 4. When the specihc values of 
the input parameters are chosen to be A = 0.01, //mp = 0.002, and A/// = 0.002, then 



















the output parameters are: a/l-p = 2.365 x 10^, n = 3.344 x 10^, Ri/lp = 1.76104695 x lO’^ 
and 6 (j){Ri)/mp = 10“®. (Here Ip is the Planck length.) The held 0(r) approaches /i for 
r > i? 2 - In the numerical integration 50(i?i) is kept hxed while Ri is varied. Fine-tuning to 
the ninth digit is necessary! If Ri is taken to be slightly bigger then 0(r) starts to deviate 
from fi in the negative direction, diverging to —cxo as r increases. If Ri is taken to be 
slightly smaller then 0 (r) starts to deviate from /i in the positive direction, heading for /2 
as r increases. With just the right value the space-time becomes nearly de Sitter outside 
the shell. The value of 6(j) at the origin (r = 0) is found from (|2.18|) to be 1.40 x 10“^^^°, 
which explains why one cannot numerically integrate 0 (r) starting from r = 0 . 

The behavior of the solution in region III is easily inferred. From the numerical inte¬ 
gration in region II both H 2 = i 7 (i? 2 ) and p 2 = p{R 2 ) are determined. In region III the 
metric can be written in the form 


H{r) = 1 - 
p{r) = p 2 . 


2GM 


(2.23) 


Here M is the mass ascribable to the shell. Once Too vanishes H{r) must take this form. 
The value of M may be determined numerically by htting Eq. (p.23|) just outside the shell. 
The location of the horizon, th, is determined by H{rH) = 0. The field 0(r) remains very 
close to /i in region HI. 

At this point due caution is required to continue the solution because the static coor¬ 


dinates dehned in (|2.2|) do not cover the whole of space-time. The global structure of the 
space-time and of the solutions is worked out in Section 4 where different coordinates are 
employed. 

The solutions we found have shell structure at a radius R whose width w is very small. 
The radius R is smaller than the horizon radius th-i but is of the same order as a. In the 
next section we shall see that R becomes smaller as / becomes smaller, but remains of 
order a even in the / ^ 0 limit. This implies that the shell structure is cosmic in size. 

An estimate of the order of magnitude of w as well as the conditions necessary for the 
existence of the cosmic shell are obtained by a simple argument. Return to Eq. (|2.11 ). 


Ln , ^ ^ AnGrHct)'^ = W[ 


Hct>"+{- 




r 


(2.24) 
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In the shell region (region II, r ~ i?) of typical solutions, H is of order one and drops 
sharply. In (|2.24|) the 2H/r term is negligible compared with the H' term. The remaining 
terms, Hcf)”, H'cj)', and are all of the same order of magnitude so that j/nP' ^ Xf^ or 


w ~ 


1 

7X7 ■ 


(2.25) 


Hence the thickness of the shell is determined by the parameters in the scalar held potential. 
The radius R is smaller than but of the same order as a. The solutions exist only if H', 
which is negative, dominates over In other words w~^ > 4:7rGR{f /w)"^. Making 

use of i? ~ a, = 9/(167rGA/^A/) and (|2.25|) , one hnds 

4aA 


47r 


mpj 


w 

< - rsj 

a 


mp 


< 1 


/ 



(2.26) 


We shall see in the next section that these relations are satished in the solutions obtained 
numerically. 


3. Numerical Analysis of the Nonlinear Regime 

In the last section we solved the linearized held equation for ^(r) in regions I and III 
where the deviation from fi is small, and we sketched the behavior in the nonlinear region 
II where the shell structure appears. In this section we present numerical results for region 
II. As discussed in the last section the boundary between regions I and II, located at the 
matching radius Ri, is rather arbitrary, subject only to the condition that the linearization 
is accurate up to that radius. Precise tuning is necessary for the pair Ri and 6(j){Ri) to 
obtain a solution to all the equations. Technically it is easier to keep 6(j){Ri) hxed and 
adjust the matching radius Ri. If Ri is chosen too small 0 comes back toward, but cannot 
reach, /i. What happens is that 0 eventually oscillates around 0 = 0 as r increases. If Ri 
is chosen too large 0 comes back to fi at hnite r to roll over it and decreases indehnitely 
toward —oo. For each sufficiently small value of 50 chosen as the matching value there 
exists a desired solution with Ri = R^. The critical value Rc can be determined numerically 
to arbitrary accuracy. 
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Figure 4: 0(r) of a solution with f/mp = 1.40 x 10 “^, A/// = 0.002, A = 0.01. 0 and 
r are in the units of mp and Ip, respectively. The maximum value of 0 is smaller than 
/2 ~ 0.0014 ■ mp. 


One example of such a solution is displayed in fig. ^for the parameters //mp = 1.40 x 
10“^, A/// = 0.002, A = 0.01 and = 1.0 x 10“^. The matching radius Ri is 

fine-tuned to ten digits: Ri/lp ~ 3.7038855228 x 10^. The shell is very thin compared to 
the radius of the shell, lying in the region 3.704 x 10^ < r/lp < 3.710 x 10^. 

In the shell region both H{r) and p(r) decrease in a two-step fashion. See hgs. and 
Inside and outside the shell H{r) is given by ( | 2 . 1 (j|) and o, whereas p(r) assumes 
the constant values 1 and 0.6256, respectively. 

The change is induced by the non-vanishing energy-momentum tensor components 
Too = —T 22 = —T 33 and Tn. They are displayed in fig. ^ and fig. ||. The energy den¬ 
sity Too has two sharp peaks associated with the rapid variation of /. The radial pressure 
Til, on the other hand, steps up quickly, remains constant within the shell, and then steps 
down again. The value of Tn is very small (~ 5 x 10“^^m/) compared with the maximum 
value of Too(~ 2 x 10“^"^m|>). The contributions of the kinetic energy, and potential 

energy, V[(()], almost cancel each other. 

In this example 0(r) makes a bounce transition /i —/2 —>■ fi once. As the numerical 
value of / becomes smaller than 0.002mp a new type of solutions with double bounces 
/i —> /2 —/i —> /2 — fi emerge. One example is displayed in £g. ^ 

The spectrum of the shell solutions depends on the parameter /. In fig. 0 R/a, 
r^/a, and R/th are plotted as functions of /. (Finding the solution numerically becomes 
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Figure 5: H{r) of a solution with f/m-p = 1.40 x 10“^, A/// = 0.002, A = 0.01. H{r) 
decreases in two steps in the shell region. H{r) in the de Sitter space (= 1 — /a^) is 

plotted in a dotted line. 

extremely difficult as / becomes smaller.) There are several features to be noted. First, 
with fixed values of A and A/// solutions exist only for / < /max. For example, with 
A = 0.01 and A/// = 0.002 we find that /max ~ 0.006. Second, as / goes to 0, rj^/a 1 
and R/a —> 0.8. The relative size R/a of the shell remains of order one even in the / —> 0 
limit. Third, as / decreases shell solutions with multiple bounces become possible. We 
cannot determine how many times 0 can bounce inside the horizon because the numerical 
evaluation becomes extremely difficult when f /rrip < 0.001. 

In the potential we are analysing, 0 = /i and 0 = /2 correspond to the false and true 
vacua, respectively. For the existence of shell solutions the fact V[fi] > V[f 2 \ is not crucial, 
however. In the example displayed in fig. ^ for instance, 0(r) swings from 0(0) fi to 
0max < f 2 , and returns to /i. Numerically V[(praax\ > td[0(O)]. As we are solving classical 
differential equations, only the form of the potential IA[0] between /i and 0max is relevant. 
The form of the potential for 0 > 0max does not matter. The second minimum can be 
higher than the first one; V[fi] < V[f 2 ] < 

Before concluding the section we would like to add that we have found no solution in 
which 0 makes a transition from /i to 02 as r varies from 0 to cxo. 
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Figure 6: p{r) of a solution with f/mp = 1.40 x 10 ^,A/// = 0.002, A = 0.01. p{r) 
decreases in two steps, from 1 to 0.6256, in the shell region. 


4. Global Structure of the Space-Time and Solutions 

In the previous section we found novel solutions to a theory with a scalar held coupled 
to gravity in static coordinates. The static coordinates in (^.2|), however, do not cover 
all of space-time. In this section we construct coordinates which allow for an extension 
of the solution to the full space-time manifold. These coordinates are smooth across the 
horizon. First we illustrate the construction with de Sitter space, and then consider the 
more general case which is applied to the shell solution. 

The X metric of the de Sitter space is given by 


ds^ = 


-dr^ + cosh^ r [dx^ + sin^ y dVP'^ 


(4.1) 


One may regard the de Sitter space as a hypersurface in the hve-dimensional Minkowski 
spacetime constrained by the condition 


2 I 2 I 2 I 2 2 I 2 

+ 2/2 + 2/3 + 2/4 — ® + 2/0 • 


(4.2) 


The metric (^7T| ) and hypersurface ([4.2| ) cover the entire de Sitter space, whereas the static 
metric covers only half of the space. Furthermore, the static metric has a coordinate 
singularity at r = a. 
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Figure 7: The energy density Too for a shell solution with f/mp = 1.40 x 10“^, A/// = 
0 .002, A = 0.01. Too and r are in the units of rrip and Ip, respectively. The energy density 
is localized in the two shells over the de Sitter background. 


The relationship among these coordinate systems are easily found. The metric (4.1) 
can be transformed to a metric conformal to the static Einstein universe. 

\—dri‘^ + dx^ + sin^ x d^'^ 


COS^ T] 


71 


^ ^ - 2tan \e , 


71 71 

- <T]< - 
2 ~ 2 


(4.3) 


Suppressing 5^, or 9 and 0 variables, one can map the whole de Sitter space to a square 
region in the y-r/ coordinates. Null geodesics are given by straight lines at 45 degree angles. 
The static and hypersurface coordinates are related by 

t 


^ q 2 _ r'i sinP . 


Vo = 


— a? cosh 


for r < a 

for r > a 


Hi = r cos 6 
y 2 = r sin 9 cos 0 
ys = r sin 9 sin 0 


t 


2/4 = 


yo? — r^ cosh — for r < a 
a 

t 


(4.4) 


— a? sinh — for r > a 
a 


14 

















Figure 8: The radial pressure Tn for a shell solution with f/m-p = 1.40 x 10“^, A/// = 
0.002, A = 0.01. Til and r are in the units of rup and /p, respectively. The pressure is 
higher and constant between the two shells. 


Similarly, ([4.2D and (|4.0|) are related by 


Hq = atan ?7 
yi = a sec rj sin y cos 6 
?/2 = a sec rj sin y sin 6 cos 0 
Us = a sec rj sin y sin 9 sin 0 

yi = a sec rj cos y. (4.5) 

The region inside the cosmological horizon (r < a, —oo < t < oo) in the static metric 
corresponds to the left quadrant in the conformal metric (0 < y < ^tt, \t]\ < 7 r /2 — y) with 
the relations 

r sin y 

a COST] 

t- = + ( 4 . 6 ) 

a 2 cos y — sin ?7 

Similarly, the region outside the cosmological horizon (r > a, —cxo < t < oo) in the static 
metric corresponds to the upper quadrant in the conformal metric (0 < y < tt, \ri\ > y—|vr) 
with the relations 

r sin y 

a cos r] 
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Figure 9: A solution with double bounces with f/mp = 1.40 x 10“^, A/// = 0.002, A = 
0.01. A solution with multiple bounces appears at a smaller shell radius than a solution 
with a single bounce. 


t 1 , sm r? + cos y 

- = - In- . 

a 2 sin ?7 — cos x 


(4.7) 


See fig. 1^. 

For the general static metric given by ( p.2|) , with H = H{r) and p = p{r), the construc¬ 
tion of the conformal coordinates analogous to (|4.3D proceeds as follows. We suppose that 
H{r) has a single zero at whereas p{r) > 0. The new radial coordinate is defined by 


r*(r) 





p{r') 

H{r') 

p{r') 

H{r') 


for r < Th 
for r > th- 


(4.8) 


It has a logarithmic singularity at th, diverging there as r* ~ —^b\n\r — th] where b = 
—2p{rH)/H'{rH) > 0. New coordinates are introduced by 


tanu = 



tann = 


(4.9) 


The upper sign is for r < rn and the lower sign is for r > th- The u-v coordinates are 


related to the U-V coordinates in the Kruskal-Szekeres (KS) coordinate system |T^ for 
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f/rn^ 


Figure 10: The /-dependence of the shell solutions. A = 0.01 and A/// = 0.002 are fixed. 
R, vh, and a are the radius of the shell, the horizon length where H{r) vanishes, and the 
horizon length in the de Sitter space, respectively. The ratios of various pairs are plotted. 
Circles, squares, and diamonds correspond to data points obtained. 


the Schwarschild solution, and to the Gibbons-Hawking (GH) coordinate system [|^ for 
de Sitter space. The connection between the KS and GH coordinate systems has already 
been discussed in [1^. Essentially tanu = U + V and tann = V — U. The metric becomes 

Ab^F{u,v) 


ds = 


cos^(m -|- v) 


dudv + r dVt 


H 

/ \ 

/ \ 


[ 1 — tanu tann 1 

(1 — cotwcott'l 

H 

/ , 2r* \ 

^ ( r < Th 


[ 1 ± cosh —— ) 

for < 

2p2 

^ b 2 

[r > Th 


F{u,v) = 

‘ip- ' ^ ^ 

. O/vi r r* r’.-T 

(4.10) 

The static metric covers the region interior to the bounding lines u = 0,u + v = ^7t, and 
u — V = ^ 71 . The horizon in the static metric, r = th, corresponds to the single point 
u = V = 0. The function F is non-vanishing and hnite there as ~ l/|r — th]- In the 
u-v coordinates the metric is regular on n = 0 so that the extension to the region u < 0 
can be made, whereas the static metric covers only half of the space. 

Applied to the de Sitter space we hnd 


th = b = a , 


^-r*/b _ 


a + r 


1/2 


F = 1 


(4,11) 
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r=oo 


71 



Figure 11: Penrose diagram of the de Sitter space. 


and the metric ([4.31) is recovered by 


TT 


r] = u + v, X = v-u + -. 


(4.12) 


The shell solution found in the preceeding sections can be extended to the entire space- 
time. In the region III (i ?2 < f') defined in the static metric H{r) is given by ( ^.23| ), and 
the location of the horizon is determined by H{rH) = 0. In region III, 5cf){r) is very small 
so that its equation of motion can be linearized. We divide region III into two; region 
Ilia (i ?2 < r < 2rH — R 2 ) and region Illb (r > 2rH — i? 2 )- In region Ilia, H{r) can be 
approximated by 

H{r) ~ /l(^l 

where 


r \ 

~ Z 2 ” / ’ 

(4.13) 

rjjJ 


:-i). 

(4.14) 


The coefficient A has been chosen such that both (|2.23|) and (|4.13|) have the same slope 
at the horizon. In the examples described in Section 3, errors caused by (|4.13|) 


are less than 15% in region Ilia. Now we write the linearized version of (|2.11|) in terms of 
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(4.15) 


Since (50(r) must be regular ai r = rn {y = 0), the solution is 

S(f){r) = (50(rH) ■ F(f + m, 1; y) 


(4.16) 


The normalization d(j){rH) must be such that d(j){r) matches at r = i ?2 with the value 
determined by numerical integration in region II. Essentially (50(r) decreases exponentially 
when R 2 < r < Th- 

Near r = R 2 the hypergeometric function behaves as 

1 


so that 


F{i + iK,f-iK,l;y) ~ 
6(j)'{r 


^-1/4(1 _ yyi/2 exp {2 k sin-1 


2r 


K 


(4.18) 


Hir) rl ^y{l-y) 

The value of the right side of ( 4.18|) at r = i ?2 can be compared with the value obtained by 
direct numerical integration in region II. In one example with f/mp = 2 x 10“^, A/// = 
2 X 10-3, A = 0.01, R 2 = 1.765 x 10^ a = 2.365 x lO^rj^ = 1.985 x 10^ k = 3344, and 
K = 4481. The numerical value for 6 (j)'{R 2 )/ 6 (j){R 2 ) is —0.000848, whereas the value from 
( [4.18| ) is —0.000985. With the uncertainty in the value of k caused by the approximation 
( [4.13| ) taken into account, one may conclude that the agreement is rather good. To gauge 
the difficulty of determining the solution numerically for all values of r, we note that 

6 (j){rH) ~ 10-1850 ■(50(i?2). 

The most important observation is that the solution 0(r) is regular at the horizon, 
r = Th, though the static metric {t, r, 6, (p) is not: r = rn in the static metric is a coordinate 
singularity. As seen above, r = th corresponds to u = n = 0 in the conformal metric ([4.10| ) 
which is non-singular. Near the horizon tan u tan v ^ r — th so that 

9 


4> = fi + S(j){rH)^l - ”^ 16 )“^"^-} ■ 


(4.19) 
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Furthermore, as tanutanu = e 2 r-.(r)/fe^ static solution (j){u,v) depends only on 
tan M tan n in the entire space-time. This implies that 

(j){—u,—v) = (j){u,v) . (4.20) 

In the Tj-x coordinates (S'^-metric) dehned by (|4.12|) 

HV, X) = tt - x) = (j){r], TT-x) ■ (4.21) 


This relation shows that there is a mirror image of the shell structure in the other hemi¬ 
sphere in the S^-metric. See £g. 12 . The shaded regions in the figure represent the shells. 
The size of the shells is invariant, as is obvious in the static coordinates. In the or con¬ 
formal metric an observer at the center (y = 0 or tt) sends a light signal toward the shell. 
Another observer at the shell, upon receiving the signal from the center, sent a signal back 
to the observer at the center. The total elapsed time, in terms of the proper time of the 
observer at the center ar where r is related to rj by ( [4.3| ) at y = 0, is a ln(a + R)/{a — R), 
independent of when the signal is emitted initially, as straightforward manipnlations show. 


71 r=oo 



Figure 12: Global structure of the shell solution in the conformal coordinates. The mirror 
image appears in the other hemisphere. 
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Finally, in the region outside the horizon in the static metric, or in the upper or lower 
quadrant in the Penrose diagram, S(j){r) decays as The proof goes as follows. For 

r 3> th, H ~ 1 —(r^/a^), so that 50(r) satishes Eq. ( 2.17 ). Hence it is a linear combination 
of Ml = F(| + m, I — in, |; z) and U 2 = + iK,j— in, z). These can be written 

as 


Ml 

M 2 


r 1 


IK, 


1-2 n 


IK 


-2 


Lr(i- 

|r(| + z/.)-^ 

F(i + 


f(| + I + ZK, 1 + 2m; 

1 - 2 m; . 


From this one can write 


0 ~ /i + H f - j sin (2k In —\- 6j 


when r ^ Th. 


(4.22) 


(4.23) 


5. False Vacuum Decay 

If the region inside the shell is in a false vacuum state one should consider its quantum 
decay to the true vacuum state. The lifetime of the false vacuum may be determined 
semi-classically using the methods of Coleman et al. without or with H gravity taken 


into account. The rate per unit volume for making a transition from the false vacuum to 
the true vacuum is expressed as 

F 


1 / 


= Ae-‘‘"'[l + 0(h)] 


(5.1) 


where Planck’s constant has been used here to emphasize the semi-classical nature of the 
tunneling rate. For the potential being used in this paper we hnd the 0(4), Euclidean 
space, bounce action (neglecting gravity) to be 

3 


B. = 


367r^ 


/ 


^ \^f) 

where the radius of the critical size bubble which nucleates the transition is 

3 [2 


(5.2) 


pc = 


a/Va- 


(5.3) 
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For any sensible estimate of the coefficient A the lifetime of the false vacuum will exceed 
the present age of the universe when the condition 


A 



< 1 


(5.4) 


is satished. This calculation is based on the thin wall approximation, which is valid when 
the critical radius is large compared to the coherence length of the potential, namely 
1/ J\V"\. This condition translates into 


6 /- 


(5.5) 


With gravity included the bounce action is 


B 


Bo 

[1 + {p,/2R)^f ' 


(5.6) 


Gravitational effects are negligible when pc < R. When this condition is not fulhlled the 
shell radius is too small to accommodate even a single nucleation bubble and therefore 
nucleation is further suppressed. 


6. Summary 


In this paper we have reported the discovery of shell-like solutions to the combined held 
equations of gravity and a scalar held with a double-minima potential. These solutions 
exist in a space that is asymptotically de Sitter. The range of parameters which allow such 
solutions are very broad. If anything like these structures exist in nature they most likely 
would have been created in the early universe and are therefore cosmological. We know of 
no other way to produce them. 

To make matters interesting, let us suppose that the cosmological constant suggested 
by recent observations of distant Type la supernovae arises from the universe actually 
being in a false vacuum state. A best ht to all cosmological data ||^ reveals that the present 
energy density of the universe has the critical value of Cc = SHq/SttG, with one-third of it 
consisting of ordinary matter and two-thirds of it contributed by the cosmological constant. 
Suppose that the cosmological constant arises from a potential of the form we have been 


22 





analyzing. Then with a present valne of the Hnbble constant of Hq = 65 km/s-Mpc we 
hnd that 


a = 


and so 


/3 1 


= 1.7 X 10^® m 


(A/^A/) = 2.4 X 10-“ eV. 


( 6 . 1 ) 


( 6 . 2 ) 


This is a constraint on the parameters of the potential, A, / and A/. Althongh they cannot 
be determined individnally from this data, we can place limits on them snch that a shell 
strnctnre might arise. Recalling (p.26|) we hnd 


/ < 30A-^/® MeV , 

w > 6.3A“^/^fm . (6.3) 


This is pnre specnlation of conrse. A cosmological constant, if it exists, may have its origins 
elsewhere. Bnt if it does arise from a false vacunm, a variety of qnestions immediately 
present themselves. Is 0 a new held, not present in the standard model of particle physics, 
whose only pnrpose is this? Where does the energy scale of 2.4 meV come from? Why 
shonld 1A[0] have a global minimnm of 0, especially when qnantnm mechanical hnctnations 
are taken into acconnt? To these qnestions we have no answers. 
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